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ABSTRACT 

With a uniform VLT SINFONI data set of nine targets, we have developed an empirical grid of J,H,K spectra of the atmospheres 
of objects estimated to have very low substellar masses of -5-20 M Sup and young ages ranging from -1-50 Myr. Most of the targets 
are companions, objects which are especially valuable for comparison with atmosphere and evolutionary models, as they present rare 
cases in which the age is accurately known from the primary. Based on the youth of the sample, all objects are expected to have low 
surface gravity, and this study investigates the critical early phases of the evolution of substellar objects. The spectra are compared 
with grids of five different theoretical atmosphere models. This analysis represents the first systematic model comparison with infrared 
spectra of young brown dwarfs. The fits to the full JHK spectra of each object result in a range of best fit effective temperatures of 
± 150-300K whether or not the full model grid or a subset restricted to lower log(g) values is used. This effective temperature range is 
significantly larger than the uncertainty typically assigned when using a single model grid. Fits to a single wavelength band can vary 
by up to 1000K using the different model grids. Since the overall shape of these spectra is governed more by the temperature than 
surface gravity, unconstrained model fits did not find matches with low surface gravity or a trend in log(g) with age. This suggests that 
empirical comparison with spectra of unambiguously young objects targets (such as those presented here) may be the most reliable 
method to search for indications of low surface gravity and youth. Based on comparison with previous observations, the SINFONI 
spectra represent a second epoch for the targets 2M0141 and DH Tau B, and the combined data show no variations in the spectral 
morphology over time. The analysis of two other targets, AB Pic B and CT Cha B, suggests that these objects may have lower 
temperatures, and consequently lower masses, than previously estimated. 

Key words, planetary systems, stars:brown dwarfs, stars:atmospheres, binariesxlose, techniques:high angular resolution 



1. Introduction 

Since the discovery of the first substellar companions GD 165B 
(Becklin & Zuckerman 1988) and Gl 229B (Nakajima et al. 
1995), there have been many surveys for additional low mass 
companions, but a very limited number have been discovered 
(e.g., McCarthy & Zuckerman 2004; Metchev & Hillenbrand 
2006), particularly at the lowest masses (Zuckerman & Song 
2009). Substellar objects that are companions to nearby stars 
and brown dwarfs are especially valuable for comparison with 
atmosphere and evolutionary models, as they present rare cases 
in which the brown dwarf age is inferred by the accurate age 
of the primary. Recent models (e.g., Fortney et al. 2008) have 
indicated the importance initial conditions may have in deter- 
mining the flux from young objects, making confirmed young 
brown dwarfs important to gauge the effects of accretion history 
and atmospheric processes at low surface gravity. 



* based on observations obtained at the Paranal Observatory , Chile 
for ESO programs 279.C-5010(A), 080.C-0590(A), 077.C-0264(A), 
078.C-0800(B), & 078.C-0800(A) 



The effects of low surface gravity on the infrared spectra of 
brown dwarfs have been investigated with isolated members of 
star-forming regions and young clusters. In the J-band, the alkali 
metal lines of K I and Na I and the FeH absorption were found 
to have lower equivalent widths in young brown dwarfs relative 
to field objects of the same spectral type in an initial sample of 
three young objects (McGovern et al. 2004) and a larger sample 
of 23 brown dwarfs in Upper Sco (Lodieu et al. 2008). Lower 
pressures in the atmospheres of the young objects cause a lower 
amount of recombination of K II and Na II to form K I and Na 
I which, along with pressure-broadening (Rice et al. 2010) and 
other processes, impacts the line strengths relative to the field 
objects, resulting in weaker lines. An index based on the Na I 
doublet was defined to assess surface gravity in the youngest 
objects (Allers et al. 2007). The triangular shape of the H-band 
was noted in the spectra of low luminosity cr Orionis members 
(Lucas et al. 2001) and associated with a decrease in H2 col- 
lision induced absorption due to low gravity in a young field 
object (Kirkpatrick et al. 2006). Rice et al. (2011) demonstrate 
how the relative opacities from CIA-H 2 absorption and H 2 at 
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different gravities shape the spectra of young objects. The CO 
bandhead shape in the K-band has also been linked with lower 
surface gravity due to the sensitivity of the absorption to micro- 
turbulence, and has been used to differentiate dwarf and giant 
stars (Kleinmann & Hall 1986). 

Given the intrinsically low luminosities and temperatures of 
substellar objects, infrared observations provide higher signal- 
to-noise spectra than optical measurements, making this wave- 
length range important for the characterization of the physical 
properties of brown dwarfs. For example, J-band spectra taken 
at two resolutions of a sample of M6-M9 objects including 15 
brown dwarfs with ages ~1-10 Myr were compared with a grid 
of model atmospheres to infer a number of properties - effective 
temperatures, surface gravities, radial velocities and projected 
rotational velocities (Rice et al. 2010). 

This paper presents the combined flux-calibrated J,H,K spec- 
tra for a set of young (2-50 Myr), low mass (5-25 Mj up ) substel- 
lar objects and compares the observations with a suite of syn- 
thetic spectra developed from five different theoretical model at- 
mosphere simulations. The range of physical properties inferred 
from the different models is discussed, along with the discrepan- 
cies obtained from fits over portions of the J,H,K wavelength 
range rather than the full combination. The spectral shape is 
the focus of this paper, and the investigation of spectral lines 
and longer wavelengths covering the L-band are the subject of a 
forthcoming paper (King et al. 201 1, in prep.). 

2. Sample and Summary of Previous Observations 

The sample consists of 9 objects estimated to have very low 
masses of ~5-20 M ]up and young ages ranging from ~l-50 Myr, 
placing them at the boundary of the brown dwarf and planetary 
mass regimes. The sample covers three age bins, with three ~1- 
3 Myr old members of star-forming regions, three members of 
the ~8 Myr TW Hydra Association, and three targets with ages 
of ~30-50 Myr. Eight of the targets are members of binary sys- 
tems discovered with adaptive optics imaging surveys, and one 
object is an isolated young brown dwarf. The youngest objects 
may still be surrounded by material from their formation envi- 
ronment and have the possibility of accretion and disks, while 
the older objects likely present bare photospheres. All are ex- 
pected to have low surface gravity. The distinct position (dered- 
dened, when necessary) on the color-magnitude diagram of this 
young, low surface gravity sample compared to field objects is 
shown in Figure 1 ; the trend is as predicted in theoretical evolu- 
tionary models (e.g., Saumon & Marley 2008). The targets show 
redder colors than field L and T dwarfs, analogous to the imaged 
planets orbiting HR 8799. For all members of the sample, there 
are archive infrared J,H,K spectra from the VLT integral field 
spectrograph SINFONI, forming a uniform data set with which 
it is possible to construct an empirical grid of the atmospheres 
of substellar objects covering the critical early phases of their 
evolution. 



2. 1 . Star-forming region systems 

High angular resolution imaging surveys targeting members of 
the nearest star-forming regions have identified a number of very 
young (~l-2 Myr) wide orbit (>100 AU) low mass brown dwarf 
companions. The pre-Main Sequence sample considered for this 
paper are the three companions with archive SINFONI J,H,K 
spectra: GQ Lup B (Neuhituser et al. 2005), DH Tau B (Itoh 
et al. 2005), and CT Cha B (Schmidt et al. 2008). The clos- 
est companion is GQ Lup B, with a separation of 077, or ~105 
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Fig. 1. Color-magnitude diagram of Field M (grey), L (green), T (blue) 
objects, the sample targets (filled red) and other young companions (un- 
filled red), and the HR 8799 planets with measured J and K magnitudes 
(black). The field data were taken from the DwarfArchive.org and re- 
stricted to objects with parallax uncertainty < 10% and color uncertainty 
<0.5 mag. The values used for the targets are given in Tables 1 and 2. 
For the young companions CHXR 73B (Luhman et al. 2006), 1RXS 
J1609B (Lafreniere et al. 2008), GSC 06214 B (Ireland et al. 2011) 
and HD 203030B (Metchev & Hillenbrand 2006) that are plotted on 
the figure, but not included in this study, the data were taken from the 
discovery papers. 

Table 2. Photometry 



Source 



H 



K 



CT Cha 


A 


9.72 4 


0.02 


8.94 4 


0.05 


8.66 4 


.0.02 




B 


16.6 


t0.3 






14.9 


±0.3 


DH Tau 


A 


9.77 i 


0.02 


8.82 4 


0.03 


8.18 4 


:0.03 




B 


15.71 


t0.05 


14.96 


± 0.04 


14.19 


±0.02 


GQLup 


A 


8.69 4 


0.04 


7.70 4 


:0.03 


7.10 4 


:0.02 




B 


14.90 


t 0.11 






13.34 


±0.13 


2M1207 


A 


13.00 


t0.03 


12.39 


± 0.03 


11.95 


±0.03 




B 


20.0 


t0.2 


18.09 


± 0.21 


16.93 


± 0.11 


TWA 5 


Aa 


8.40 4 


0.07 


7.69i 


:0.04 


7.39 4 


:0.04 




Ab 


8.5 4 


0.2 


7.79 4 


:0.05 


7.62 4 


-0.08 




B 


12.6 


t0.2 


12.14 


± 0.06 


11.4 


± 0.2 


AB Pic 


A 


7.58 4 


0.03 


7.09 4 


0.03 


6.98i 


0.03 




B 


16.18 


±0.10 


14.69 


±0.10 


14.14 


±0.08 


GSC08047 


A 


9.06 4 


0.03 


8.53i 


:0.06 


8.41 4 


:0.03 




B 


15.9 


±0.1 


15.45 


±0.02 


14.8 


±0.1 


2M0141 




14.83 


±0.04 


13.88 


±0.03 


13.10 


±0.03 



AU at a distance of 150 ± 20 pc (Franco 2002), and has been 
the subject of considerable prior investigation into the inferred 
companion mass and evidence for ongoing accretion. Based on 
a comparison of theoretical models with photometry and K- 
band spectroscopy, the mass was initially estimated as 1-42 M Jup 
(Neuhauser et al. 2005), while subsequent analysis of a larger 
wavelength range (e.g., Marois et al. 2007) revised the mass 
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Table 1. Sample 



OUUILC 


Pri m a r v 
i Hindi y 


Comp 


A op 


M - 

lyI prim 


M 


J3 


Sep. 


lVf pmhpr^hin 

IVlClllL't;! Jill LJ 


A , 


Refs. 








(Myr) 


(M ) 


(M Jup ) 


(pc) 


(AU) 








GQ Lup B 


K7 


L1.5 


1 


0.7 


24±12 


150±20 


100 


Lupus SFR 


0.1+0.2 


1-5 


UH lau a 


MU.o 


LI 


i 


V.iS 


11± 3 


14t)± 1 J 




Taurus SFR 


0.J+U.3 


6,2 


CT Cha B 


K7 


>M8 


2 


0.7 


17±6 


165±30 


440 


Chameleon SFR 


0.4+0.2 


7-8 


2M1207 A&B 


M8 


L5 


8 


25±5 M Jup 


6-10 


52.4+1.1 


46 


TW Hydra Assoc. 




9-12 


TWA 5 B 


M1.5 


M8/8.5 


8 


0.40 


-20 


44±4 


98 


TW Hydra Assoc. 




13-15 


AB Pic B 


K2 


L0-L1 


30 


0.84 


13-14 


47.3+1.8 


248 


Tuc-Hor Assoc. 




16-17 


GSC08047 B 


K3 


M9.5 


30 


0.80 


25+10 


-85 


279 


Tuc-Hor Assoc. 




18-19 


2M0141 


LO 




1-50 


15±10M Jup 




-35 




Tuc-Hor//? Pic? 




20 



References: (1) Neuhauser et al. (2005) (2) Luhman et al. (2006) (3) Marois et al. (2007) (4) Franco (2002) (5) Batalha et al. (2001) (6) Itoh et al. 
(2005) (7) Schmidt et al. (2008) (8) Luhman (2004) (9) Gizis (2002) (10) Chauvin et al. (2004) (11) Chauvin et al. (2005a) (12) Mohanty et al. 
(2007) (13) Lowrance et al. (1999) (14) Webb et al. (1999) (15) Konopacky et al. (2007) (16) Chauvin et al. (2005c) (17) Bonnefoy et al. (2010) 
(18) Chauvin et al. (2005b) (19) Neuhauser & Guenther (2004) (20) Kirkpatrick et al. (2006) 



range to 10-20 Mj up . In Table 1, the mass is listed as 24 ± 12 
M j up (Luhman et al. 2006), since this estimate applied the same 
approach to both GQ Lup B and another target DH Tau B. The 
SINFONI spectra reported in Seifahrt et al. (2007) suggested ev- 
idence for active accretion from the Pa /3 emission line, while 
independent J,H,K spectra from OSIRIS (McElwain et al. 2007) 
and NIFS (Lavigne et al. 2009) do not confirm the line, which 
may have resulted from contamination from the primary. The 
SINFONI spectrum of GQ Lup is included in this paper for com- 
pleteness. 

The companion DH Tau B (Itoh et al. 2005) is in a wider 
orbit, at 273 (330 AU) from the primary, making contamination 
from the host star unlikely. The estimated mass for the compan- 
ion, 1 1^1° Mj up , listed in Table 1 is based on the same model 
comparison applied to GQ Lup B (Luhman et al. 2006) and 
within the 5-40 M Sup range reported in Itoh et al. (2008). The pre- 
vious CISCO J,H,K spectrum (R-440) showed absorption from 
H 2 bands indicating a low effective temperature, and no evi- 
dence of emission lines from accretion was present (Itoh et al. 
2005). The SINFONI data presented in this paper were taken at 
higher resolution (R-1500) than the CISCO spectrum and at a 
different epoch, enabling a search for variability. 

The widest young companion CT Cha B is 277 (440 AU) 
from the primary (Schmidt et al. 2008). Based on fitting the 
SINFONI J,H,K spectra with Drift-PHOENIX models to deter- 
mine the temperature and photometry to estimate the luminosity, 
the companion mass was estimated to be 17 ± 6Mj up (Schmidt 
et al. 2008) from comparison with evolutionary models (Baraffe 
et al. 2003; Chabrier et al. 2000; Burrows et al. 1997). The spec- 
trum shows Pa j3 emission, an indication of ongoing accretion, 
though another accretion diagnostic Br y is not seen (Schmidt 
et al. 2008). In this paper, we reconsider the SINFONI spectra 
in comparison with a larger set of models and with a different 
treatment of extinction. 

2.2. 71/1/ Hydra Association systems 

Two members of the TW Hydra Association - TWA 5 and 
2M1207 - include three substellar objects, with both late-M and 
mid-L spectral types. TWA 5B was one of the earliest brown 
dwarf companions identified (Lowrance et al. 1999; Webb et al. 
1999) and is located 270 (-100 AU) from a close pair of early 
M-stars with a -5 yr orbit (Konopacky et al. 2007). Although 
the photometry of TWA 5B has been explored (Neuhauser et al. 



2000), the H- and K-band spectra have been published sepa- 
rately (Neuhauser et al. 2009). The second TW Hydra mem- 
ber, 2M1207, is composed of two substellar objects (Chauvin 
et al. 2004, 2005a) - the primary is an M8 brown dwarf (Gizis 
2002) similar to TWA 5B and the secondary is a L-type plan- 
etary mass object (Chauvin et al. 2004; Mohanty et al. 2007). 
The complete J,H,K spectrum of 2M1207B has been published 
(Patience et al. 2010), but not the corresponding spectrum of the 
primary. When compared with DUSTY models, the companion 
spectrum presents a spectral shape consistent with temperatures 
higher than expected from the luminosity (Mohanty et al. 2007; 
Patience et al. 2010). The discrepancy between spectral shape 
and flux level was initially explained by obscuration from a grey 
disk (Mohanty et al. 2007), and then updated to include mod- 
els incorporating an improved treatment of clouds (Skemer et al. 
2011) and both clouds and non-equilibrium chemistry (Barman 
et al. 2011) to simultaneously fit the spectrum shape and flux 
level. The 2M1207 B spectrum represents an example of how 
measurements of young, low mass object atmospheres can ad- 
vance the understanding of the atmospheric physics. 

2.3. Tucana-Horolgium systems 

An additional two substellar companions - AB Pic B and 
GSC 080647 B - are members of the older -30 Myr Tucana- 
Horologium Association. AB Pic B is a wide orbit 575 (-260 
AU) common proper motion companion (Chauvin et al. 2005c) 
to a K2 primary (Perryman et al. 1997). Based on a comparison 
of the system photometry with evolutionary models (Burrows 
et al. 1997; Baraffe et al. 2002) at 30 Myr, the mass of the com- 
panion is estimated to be 13-14 M Jup (Chauvin et al. 2005c), 
placing AB Pic B at the demarcation of planetary masses and 
brown dwarf masses. The J,H,K SINFONI spectra have been 
used to estimate a spectral type for the companion of L0-L1 
and to compare each band of the spectrum with theoretical at- 
mosphere models (Bonnefoy et al. 2010). In this paper, a model 
comparison of the combined flux-calibrated J-K spectrum is pre- 
sented, and the combined spectrum shows differences with the 
models not apparent in the single band comparisons. Another 
Tucana-Horologium pair with very similar spectral types to AB 
Pic is GSC 08047, which is composed of a K3 primary with a 
confirmed substellar companion -373 (279 AU) from the host 
star (Chauvin et al. 2005b). Previous spectra in the H- and 
K-bands taken with VLT/ISAAC and VLT/NACO were com- 
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Table 3. SINFONI Observations 



Source 


Grating 


Scale 


NDIT 


DIT 


#exp. 


Dates 


ESO Project ID 


PI 


CT Cha B 


J 


lOOmas 


2 


138.5s 


6 


18May07 


279.C-5010(A) 


Schmidt 


CT Cha B 


H+K 


lOOmas 


2 


138.5s 


6 


16May07 


279.C-5010(A) 


Schmidt 


DH Tau B 


J 


25mas 


1 


300s 


32 


18Dec07 


080.C-0590(A) 


Rojo 


DH Tau B 


H+K 


25mas 


1 


300s 


8 


22Oct07 


080.C-0590(A) 


Rojo 


GQ Lup B 


J 


25mas 


1 


300s 


27 


18Sep06 


077.C-0264(A) 


Neuhauser 


GQ Lup B 


H 


25mas 


1 


300s 


11 


24Ap06 


077.C-0264(A) 


Neuhauser 


GQ Lup B 


K 


25mas 


1 


300s 


8 


17Sep05 


077.C-0264(A) 


Neuhauser 


TWA 5B 


J 


25mas 


1 


120s 


6 


12Dec07 


080.C-0590(A) 


Rojo 


TWA 5B 


H+K 


25mas 


2 


60s 


7 


13Dec07 


080.C-0590(A) 


Rojo 


2M1207A+B 


J 


lOOmas 


1 


600s 


22 


25Feb07 


078.C-0800(B) 


Thatte 


2M1207A+B 


H+K 


lOOmas 


1 


300s 


24 


28Jan07, 07Feb07 


078.C-0800(B) 


Thatte 


AB Pic B 


J 


25mas 


1 


300s 


40 


5-7Dec07 


080.C-0590(A) 


Rojo 


AB Pic B 


H+K 


25mas 


1 


300s 


8 


1 1Nov07 


080.C-0590(A) 


Rojo 


GSC 08047 B 


J 


25mas 


1 


300s 


24 


6&18Jan08 


080.C-0590(A) 


Rojo 


GSC 08047 B 


H+K 


25mas 


1 


300s 


27 


5&9Jan08 


080.C-0590(A) 


Rojo 


2M0141 


J 


250 mas 


1 


300 


4 


19Oct06 


078.C-0800(A) 


Thatte 


2M0141 


H+K 


250 mas 


1 


60s 


4 


19Oct06 


078.C-0800(A) 


Thatte 



pared with late-M objects to estimate a spectral type of M8 ± 
2 (Neuhauser & Guenther 2004) and M9.5 ± 1 (Chauvin et al. 
2005b), slightly earlier than the AP Pic companion. The new 
SINFONI spectra include the J-band and the IFU observations 
are not subject to wavelength-dependent effects of differences in 
the slit centering for target and standard that affected some of the 
previous spectra. 

2.4. Young Isolated Field Object 

As part of the large scale search and characterization program for 
L-dwarfs with 2MASS, an isolated brown dwarf with indications 
of youth was identified with spectroscopy spanning 0.6 - 2.5 fim 
(Kirkpatrick et al. 2006). In the infrared, there is a J-K SpeX 
spectrum with similar resolution (R~ 1200) to the SINFONI data, 
and a NIRSPEC J-band spectrum with higher (R-2500 ) resolu- 
tion, making 2M0141 an ideal case to test for variability such as 
has been reported for TWA 30 (Looper et al. 2010). For 2M0141, 
the spectral features of weak alkali lines such as K I in the J- 
band and the triangular shape of the H-band caused by a de- 
crease in H2 collision induced absorption provided evidence of a 
low surface gravity object (Kirkpatrick et al. 2006). The effective 
temperature and surface gravity were derived from fits to atmo- 
spheric models generated with the PHOENIX code (Kirkpatrick 
et al. 2006), and r eff and log(g) were compared with evolutionary 
models (Baraffe et al. 2002) to infer an age range of 1-50 Myr 
and object mass of 6-25 Myr (Kirkpatrick et al. 2006). Although 
2M0141 currently has a less certain age, future measurements 
of the parallax and space motion may determine membership in 
the Tucana-Horologium Association or the ft Pic moving group. 
In this paper, the SINFONI spectra are compared with a larger 
set of atmosphere models and with the earlier SpeX spectrum to 
search for variability. 

3. SINFONI Observations 

The observations were taken with the SINFONI (Spectrograph 
for INtegral Field Observations in the Near Infrared) instru- 
ment (Bonnet et al. 2004; Eisenhauer et al. 2003; Thatte et al. 
1998), an AO-equipped integral field spectrograph mounted at 
the Cassegrain focus of the VLT (UT4). The target low mass 



companions are separated by less than a few arcseconds from 
the primaries, making the primaries ideal sources for the AO cor- 
rection. The observing sequence for all targets included a set of 
exposures on the target with offset positions on the array. For 
some targets, separate blank sky exposures were also recorded. 
Standard calibration observations of an early-type or solar-type 
star followed each target sequence to correct for telluric features 
and the instrument response. Additional calibrations were ob- 
tained to measure the wavelength scale, the detector dark cur- 
rent, distortion, and quantum efficiency variations. The data for 
this project were drawn from the ESO archive, and a summary 
of the instrument configuration, observation details and project 
ID codes for each target is given in Table 3. In some cases, only 
the subset of higher quality data were used; Table 3 lists only the 
observations used in this analysis. 

There are four gratings available in SINFONI - J, H, K, 
and H+K - which produce spectral resolutions of approximately 
2000, 3000, 4000, and 1500, depending on the scale. For all tar- 
gets except GQ Lup B, a combination of J and H+K observations 
were taken. For GQ Lup B, separate J, H, and K spectra were 
recorded. Among the three available spaxel (spatial pixel) scales 
of 25mas, lOOmas, and 250mas, the lOOmas scale was employed 
for the CT Cha B and 2M1207A+B observations. The lOOmas 
scale provides a field-of-view of 3'.'0 x 3'.'0, and both the brown 
dwarf primary and planetary mass secondary for 2M1207 were 
included within each observation. For the remaining target, the 
spaxel scale of 25 mas was used and the field-of-view at this scale 
is 078 x 078 which excluded the primary host star from the ob- 
servations. For one target, GQ Lup B, the diffraction spike from 
the primary intersected the position of the target companion in a 
subset of the observations. 



4. Data Analysis 

4.1. SINFONI data reduction 

The initial steps of the data reduction for all targets were per- 
formed with the Gasgano (Gebbinck et al. 2007) implementation 
of the SINFONI data reduction pipeline (Modigliani et al. 2009; 
Dumas et al. 2007). The processing began with the raw data and 
calibration files in the ESO archive and is independent of previ- 
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Fig. 2. The evolution in surface gravity over the 1-100 Myr age range 
for objects with masses of ~ 1-50 Mj up . The values are taken from the 
DUSTY00 models (Chabrier et al. 2000) and show the systematic in- 
crease in log(g) over time as the objects contract. These tracks are used 
to determine the restricted range of surface gravity to allow in model 
fits to our young targets. 



ously reported results summarised in Section 2. This approach 
was taken to ensure that all data used the same software ver- 
sion and were reduced in a consistent way. Corrections for the 
dark current, linearity response, and bad pixels were measured 
and applied to all data files. The instrument optical distortion 
was measured with a set of fibers illuminating positions across 
the field and the non-linear wavelength dispersion was measured 
with observations of an arc lamp; both the spatial and spectral 
calibrations were applied to the data with Gasgano routines. The 
calibrated two-dimensional data for each observation was con- 
verted into a three-dimensional data cube consisting of an image 
at each wavelength slice. Across the entire wavelength range, 
a stack of ~700-2000 wavelength slices images were recorded, 
however the faintest target 2M1207B required binning in the J- 
band to produce a spectrum with a higher signal-to-noise ratio. 

After corrections for the detector and optics, the sky back- 
ground emission was subtracted from the offset position cover- 
ing either blank sky or including the target at a different position 
on the array in the ABBA pattern of observations. The individual 
sky-subtracted cubes were aligned based on centroiding the tar- 
get in the top wavelength slice and then combined. To ensure the 
final spectrum traced a straight path through the data cube, each 
wavelength slice of the final data cube was shifted and aligned, 
and then the final spectrum was extracted. For the one target with 
both components on the array, 2M1207, it was necessary to sub- 
tract the halo of the primary before extracting the spectrum of the 
companion. Further details of the subtraction have been reported 
previously (Patience et al. 2010). 

The standard stars observed after the targets were used to 
correct for telluric absorption and the wavelength response of 
the detector by dividing each target spectrum by its associated 
standard after interpolating the standard spectrum over intrin- 
sic features. The shape of the standard star continuum was re- 



Fig. 3. The extinction Aj of the secondary as a function of the Aj 
of the primary for spatially resolved binary systems in Lupus (red), 
Chameleon (light blue), and Taurus (dark blue) based on observations 
of the individual components of each system (Brandner & Zinnecker 
1997; White & Ghez 2001; Prato et al. 2003). These binaries have com- 
parable separations to the target systems and are located in the same 
star-forming regions. The data show a trend of comparable extinction 
for each component. 



moved by multiplication of a black body with the same effec- 
tive temperature as the standard. To flux calibrate the spectra, it 
was necessary to rely on photometric measurements presented 
in Table 2 which were not contemporaneous with the spectra. 
The observed spectra and the spectrum of Vega (Mountain et al. 
1985; Hayes 1985) were convolved with the 2MASS response 
curves (Carpenter 2001) to calculate the scaling factor required 
to match the reported magnitude in a given filter. For the H+K 
spectra, the H-band photometry was used for flux calibration, 
because the SINFONI data covers the full 2MASS H-band filter 
bandpass. For one target, GQ Lup B, there is no reported H-band 
magnitude and the H and K spectra were observed separately. In 
this case, the H-band magnitude was estimated by assuming the 
same H-K color as DH Tau B, since GQ Lup B and DH Tau 
B both have a very similar J-K color and similar young ages. 
Another target, CT Cha B, also had no H-band measurement, 
and, in this case, the K-band magnitude was used to flux cali- 
brate the H+K spectrum. 

4.2. Comparison with theoretical model grids 

The target spectra were compared with five sets of theoreti- 
cal models, spanning a range of effective temperature and sur- 
face gravity, as described in Section 5. The model spectra were 
smoothed to the same resolution as the data and interpolated to 
the same dispersion before performing a fitting procedure using 
the least-squares statistic defined in Mohanty et al. (2007) and 
allowing the model flux to be scaled to find the best fit. This 
flux scaling is related to the object radius, enabling an estimate 
of the radius in addition to the effective temperature and sur- 
face gravity. Since each wavelength step in the SINFONI spec- 
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Fig. 4. For a fixed value of log(g)=4.5, the solar metallicity model spectra at temperatures ranging from 1400K (upper left) to 2900K (lower right) 
in steps of 300K, corresponding to the target effective temperatures. From top to bottom in each panel, the models are: DUSTY (grey), BT-Settl 
(green), Drift-PHOENIX (red), Gaia-Dusty (purple), and Marley et al. (blue). The Marley et al. models do not include temperatures above 2400K, 
and the plots show models with f se d=l, the dustiest set of models within this model grid. 

Table 4. Cloud Models" 



Physics 



DUSTY 



BT-Settl 



Drift-PHOENIX 



Gaia-Dusty 



Marley et al. 



30 dust species 
S = 1 
none 



updated DUSTY 
S = 1 
none 



none 

„-3.5 



f(a) ec a~ i5 , (ISM) 
0.00625-0.24 yum 
equilibrium 



none 

„-3.5 



5 dust species 
5 = 1 
diffusive'' 
parameterised in f sec i 
log-normal f(z, s) e 

IQ- 5 - 1 fim 
equilibrium except 

CO/CH 4 , N 2 /NH 3 

Notes: a PHOENIX (Hauschildt et al. 1999) DUSTY (Allard et al. 2001), BT-Settl (Allard et al. 2003, 2010), Drift-PHOENIX (Helling et al. 
2008; Witte et al. 2009), Gaia-Dusty (Rice et al. 2010), Marley et al. (Ackerman & Marley 2001; Marley et al. 2003; Saumon & Marley 2008; 
Freedman et al. 2008) * S quantifies the supersaturation ratio. 5 = 1 for thermal equilibrium between the gas and dust. c Parameterised in k zz . c 

f(a, z) is the grain size distribution at height z and provides the number of grains of size a. e f(z, s) is the grain size distribution for a grain species 

s at height z- ' Updated chemistry treatment compared to DUSTY. 



Opacity source 
Gas-dust phase* 
Mixing 

Settling/Rainout 
Dust sizes 

Chemistry 



43 dust species 
S = 1.001 
timescale from RHD 
timescale comparison 
timescale comparison 
0.00625-0.24 /an 
equilibrium 



7 dust species 
S =f(z,s) 
damping of turnover timescale 
equation of motion 
f(a,z) d 
10~ 6 - 100 /an 
equilibrium^ 



f(a) ec (ISM) 
0.00625-0.24 yum 
equilibrium^ 



tra is nearly constant over the JHK range, we do not employ the 
weighting technique of Cushing et al. (2008) which was devel- 
oped for spectra covering a much larger range of wavelengths. 

Two sets of fits were performed to the model grids - a fit 
including the full range of log(g) values and a fit restricted to 
log(g)>3.0 and with an upper limit based on evolutionary models 
(Chabrier et al. 2000) shown in Figure 2. The targets have known 
young ages, so the predicted evolution in surface gravity was 
used to set log(g) limits: 4.0 for the companions to pre-Main 
Sequence stars, 4.5 for the TW Hydra members, and 5.0 for the 
30-50 Myr targets. 

A final factor affecting the slope of the observed spectrum 
is the amount of extinction to each target. This is particularly 
important for the three objects in star-forming regions. None of 
the secondaries have measured values of extinction, however, the 
primaries have reported Ay values. Rather than treat extinction 
as a free parameter, we use the results of binary star component 
extinction comparisons to gauge the likelihood of very different 
Aj values for the primary and secondary. The secondary Ay is 



shown as a function of primary Aj in Figure 3 for binary stars 
in the same regions as the targets (Brandner & Zinnecker 1997; 
White & Ghez 2001; Prato et al. 2003), and a clear trend of sim- 
ilar extinction values is seen. Based on the binary stars, the tar- 
gets in this sample are assumed to have the same extinction as 
the primary, which is reported in Table 1 . 



5. Theoretical Models 

The data were compared with five grids of theoretical atmo- 
sphere models. Four of the models represent different extensions 
of the PHOENIX (Hauschildt et al. 1999) code and incorporate a 
range of dust properties and physical effects. In chronological or- 
der of development, the PHOENIX-based models considered in 
this paper include: DUSTY (Allard et al. 2001), BT-Settl (Allard 
et al. 2003, 2010), Drift-PHOENIX (Helling et al. 2008; Witte 
et al. 2009), and Gaia-Dusty (Rice et al. 2010), with the latter 
three representing ongoing modelling efforts. The earlier genera- 
tion DUSTY model is included for comparison, although the BT- 
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Table 5. Model Grid Parameters 



Model 


Teg 




log(g) 


AA 






range 


step 


range 


step 






DUSTY 


500-4000 


100 


3.5-6.0 


0.5 


2k 


2.5 - 


BT-Settl" 


400-3000 


100 


3.5-5.5 


0.5 


0.2 A 




Drift-PHOENIX 


1000-3000 


100 


3.0-5.5 


0.5 


2 A 


1.5 - 


Gaia-Dusty 


1400-3500 


50 


3.0-6.5 


0.1 


0.5 A 




Marley et al* 


900-2400 


100 


4.0-5.0 


0.25 


R~50,000 


0.5 - 


Note: " for 7/ eff 


=2000-3000K, the log(g) range is 0.0-5.5 


* the 




wa 


velength step A/1 is a function of A 




- 



1.2 1.4 



Wavelength ((j,m) 



2 2.2 24 



Settl and Gaia-Dusty models have superseded the DUSTY mod- 
els with improved treatments of abundances and line lists. The 
independent model grid developed by Marley et al. (Ackerman 
& Marley 2001; Marley et al. 2003; Saumon & Marley 2008; 
Freedman et al. 2008) is also compared with the target spectra. 

A summary of the physics incorporated into the five mod- 
els is given in Table 4; details of the models are given in the 
references listed above and the comparison of models with two 
benchmark test cases is reported in (Helling et al. 2008). The 
treatment of some physical processes such as convection is iden- 
tical for all models, while the approach to other factors such as 
condensate settling and chemistry varies between the models. 
Since dust plays a large role in determining the temperature- 
pressure structure of the atmospheres in the temperature regime 
appropriate for the targets, the number of dust species and dust 
size distributions involved in the calculations are listed, along 
with notes on the mixing and settling included in the models. 

An example of the synthetic spectra produced by each model 
over the same range of for a fixed value of log(g) = 4.5 is 
given in Figure 4. This value of log(g) was selected because it is 
within the range predicted by evolutionary models for the ages 
of the targets. For the Marley et al. spectra, the value of f sel j is 
1 which represents the dustiest atmospheres. At the coolest tem- 
perature shown, 1400 K, there is a significant difference in the 
synthetic shape of the H-band spectrum and the overall slope for 
the different models. These effects are reduced by 1700K, but 
still impact the peak of the J-band. By a temperature of 2000 K, 
the slope across the J-band ranges from flat to increasing. At hot- 
ter temperatures, the line strengths differentiate the models more 
than the overall shape, and the atomic and metal hydride lines 
will be investigated for this sample in a future paper (King et al. 
2011, in prep.). 

These model differences are primarily caused by the differ- 
ent treatments of condensate settling and the differences in dust 
opacity. Typically, a higher overall dust opacity results in a rel- 
atively lowered J-band flux and weakens the water bands while 
smoothing their edges. The H band shape is dominated by the 
two adjacent water bands and the collision induced absorption 
(CIA) which both grow stronger with decreasing dust opacities 
as flux emerges from the deeper, hotter atmosphere resulting in 
a bluer spectrum. 

Another difference between the models is apparent upon 
close comparison of the strong absorption lines and bands be- 
tween models calculated at different resolutions. In regions of 
strong, varied absorption, models calculated at lower resolution 
(2 A) look 'noisier', i.e. a larger spread in the flux values over 
small wavelength ranges (see Fig. 4). For example, the DUSTY 
models show the largest 'noise' of the spectrum over the full 
temperature range, while the Marley et al., BT-Settl, and Gaia- 
Dusty models are typically much narrower in the range of flux 
at a given wavelength. These differences are a numerical artifact 




Wavelength (\im) 



Wavelength (\im) 

Fig. 5. SINFONI J,H,K spectra of the sample, separated into three age 
groups: (top panel, top to bottom) the 1-3 Myr objects DH Tau B, GQ 
Lup B, and CT Cha B; (middle panel, top to bottom) the 8 Myr objects 
2M1207A, TWA 5B and 2M1207B; (bottom panel, top to bottom) the 
30-50 Myr objects GSC 08047B, 2M0141, and AB Pic B. For each age 
group, the spectra are ordered from hottest (top) to coolest (bottom), and 
an offset has been added to the flux for ease of comparison. The spectra 
of the youngest subset have been dereddened based on the extinction 
measured for the primaries. 



induced when the convolution with a Gaussian instrument dis- 
persion profile is undersampled, resulting in lower fidelity model 
spectra for DUSTY and Drift-PHOENIX. The numerical origin 
of the difference was confirmed by calculating one grid (Gaia- 
Dusty) with two initial resolutions and comparing the results. 
Due to the fact that the model fits in this paper are based only 
on the spectral shape, the additional width in the DUSTY and 
Drift-PHOENIX models does not impact this analysis, however 
the future comparison of spectral lines and indices requires mod- 
els with sufficient resolution to fit line strengths. 

Table 5 shows the range of T eS and log(g) covered by the so- 
lar metallicity models. The Marley et al. grids also incorporate 
a parameter f se d describing the sedimentation efficiency, and the 
values considered for f sed in this grid are 1 and 2, with the higher 
value associated with more efficient sedimentation from more 
rapid particle growth and larger average particle sizes. In addi- 
tion to the solar metallicity models summarized in Table 5, the 
Drift-PHOENIX grid also covered non-solar metallicities with 
values of -0.6, -0.3, 0.0, and +0.3. The data are first compared 
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Table 6. Best Fitting Atmopshere Model Parameters, Full Grid Comparison 



Wavelength Range 


DUSTY 


BT-Settl Drift-PHOENIX 


Gaia-Dusty 


Marley et al. 






7-eff (K) 


logg 


T eS (K) 


log# T eB (K) 




7W (K) 


logg 


T eS (K) 


logg 


fsed 


1-3 Myr Objects 


DH Tau B" 


JHK(1.10-2.46/an) 


2600 


4.0 


2400 


2.5 2000 


4.0 


2350 


3.7 


>2400* 


4.0 


2 


J(1.10-1.34/im) 


2400 


4.0 


2200 


3.5 1600 


3.5 


2350 


4.1 


>2400* 


4.0 


2 


HK (1.46-2.46//m) 


2600 


3.5 


2500 


3.5 1900 


5.5 


2550 


3.8 


>2400* 


4.0 


2 


GQ Lup B" 


JHK(1.10-2.46jum) 


1900 


6.0 


1800 


5.5 1700 


3.5 


2050 


5.5 


1600 


4.5 


1 


J(1.10-1.34jum) 


2500 


4.0 


2400 


3.5 2600 


4.0 


2500 


4.8 


>2400* 


4.0 


2 


HK(1.46-2.46yum) 


1800 


6.0 


1700 


4.5 1700 


5.0 


1800 


5.5 


1500 


4.5 


1 


CT Cha B" 


JHK(1.10-2.46//m) 


1800 


4.5 


1600 


4.0 1700 


3.5 


1950 


4.1 


1500 


4.5 


2 


J(1.10-1.34yum) 


1900 


3.5 


1400 


5.0 1700 


3.0 


2000 


3.9 


1100 


5.0 


2 


HK (1.46-2.46yCim) 


1900 


6.0 


2200 


1.0 1800 


5.0 


1950 


5.5 


1700 


4.5 


1 


8 Myr Objects 


2MASS 1207 A 


JHK(1.10-2.46/mi) 


3100 


5.5 


2800 


5.0 >3000" 


5.0 


2900 


5.3 


>2400* 


5.0 


1 


J(1.10-1.34yum) 


2900 


5.0 


2700 


4.5 >3000* 


5.0 


2850 


5.1 


>2400* 


4.0 


2 


HK (1.46-2.46yum) 


3100 


5.5 


2800 


5.5 >3000 6 


5.5 


2950 


5.5 


>2400* 


5.0 


1 


TWA 5 B 


JHK(1.10-2.46/mi) 


2500 


3.5 


2400 


2.5 2200 


3.0 


2450 


3.3 


>2400* 


4.0 


2 


J(1.10-1.34^m) 


2700 


4.5 


2600 


4.0 2800 


4.5 


2700 


5.0 


>2400* 


4.0 


2 


HK (1.46-2.46ytim) 


1800 


6.0 


2200 


1.0 1800 


5.5 


1900 


5.5 


1600 


4.5 


1 


2MASS 1207 B 


JHK(1.10-2.46/mi) 


1600 


3.5 


1500 


3.5 1500 


5.0 


1650 


4.3 


1100 


5.0 


1 


J(1.10-1.34yum) 


1700 


4.5 


2000 


0.5 1700 


4.0 


1950 


4.0 


1500 


4.0 


1 


HK (1.46-2.46yum) 


1600 


4.5 


1500 


3.5 1500 


3.0 


1650 


4.7 


1300 


4.0 


1 


30-50 Myr Objects 


GSC 08047 B 


JHK(1.10-2.46/mi) 


2200 


3.5 


2200 


1.0 1900 


4.5 


2250 


3.1 


2300 


4.0 


2 


J(1.10-1.34yum) 


2600 


5.0 


2200 


3.5 1600 


3.5 


2300 


4.2 


>2400* 


4.0 


2 


HK (1.46-2.46ytim) 


1800 


6.0 


1700 


4.5 1700 


4.0 


1850 


5.5 


1600 


4.5 


1 


2MASS 0141 


JHK(1.10-2.46//m) 


1800 


5.0 


2000 


3.0 1700 


4.0 


2000 


5.5 


1600 


4.5 


1 


J(1.10-1.34yum) 


2300 


4.0 


2200 


3.5 1600 


3.5 


2300 


4.9 


>2400* 


4.0 


2 


HK (1.46-2.46yum) 


1800 


6.0 


2100 


1.0 1800 


5.0 


1900 


5.5 


1600 


4.5 


1 


AB Pic B 


JHK(1.10-2.46/mi) 


1800 


6.0 


1600 


3.5 1600 


4.5 


1850 


5.5 


1400 


4.5 


1 


J(1.10-1.34yum) 


2400 


4.5 


2300 


4.0 1600 


3.5 


2400 


4.9 


>2400* 


4.0 


2 


HK (1.46-2.46ytim) 


1900 


6.0 


2300 


1.5 1800 


5.0 


1950 


5.5 


1700 


5.0 


1 



Note: " These star formation region members have been dereddened as described in Sect. 4. 6 Note that this T cft falls at the edge of the grid used. 



with the full grid summarized in Table 5 and then with a re- 
stricted subset with limits on the log(g) values appropriate to 
the ages of the targets, based on evolutionary models (Chabrier 
et al. 2000). Finally, the impact of metallicity is explored with 
a comparison of the data with the expanded grid with different 
metallicities for the Drift-PHOENIX model. 

6. Results and Discussion 

6.1. SINFONI Spectra 

The J,H,K flux-calibrated SINFONI spectra for the nine young 
substellar objects in the sample are plotted in Figure 5, grouped 
in the three age bands and ordered in decreasing temperature 
determined from the model comparisons described in Section 

6.2. Approximately half of the SINFONI data are previously un- 
published: all wavelength bands of DH Tau B, 2M1207B, and 
2M0141 are newly reported along with the J- and H-bands of 
TWA 5B and GSC 080647 B. For the youngest targets, the spec- 
tra have been dereddened according to the reddening law of 



Cardelli et al. (1989) using the extinction measured for the pri- 
mary, as described in Section 4. 

The overall spectral shape is dominated by the H 2 absorp- 
tion at both sides of each wavelength band and the CO absorp- 
tion band at the red edge of the K-band. All the object spectra 
exhibit a distinct triangular shape in the H-band portion of the 
spectrum, a signature of young substellar objects. As noted in 
other spectra of young objects, the peak of the K-band also oc- 
curs at a longer wavelength near 2.3/im than is typical for older 
field L-dwarfs (Mohanty et al. 2007). Two of the coolest objects, 
2M1207 B and AB Pic B, show an unusually depressed J-band 
portion of the spectrum relative to the spectra of the other sub- 
stellar objects of the same age. 

6.2. Comparison with Atmosphere and Evolutionary Models 
6.2.1. Full Model Grid 

Each object spectrum was compared with a grid of the five solar- 
metallicity models summarized in Section 5 - DUSTY, BT-Settl, 
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Table 7. Best Fitting Atmopshere Model Parameters, Restricted Grid Comparison 



Wavelength Range 


DUSTY 


BT-Settl Drift-PHOENIX 


Gaia-Dusty 


Marley et al. 






7-eff (K) 


logg 


T eS (K) 


log# T eB (K) 




7W (K) 


logg 


T eS (K) 


logg 


fsed 


1-3 Myr Objects 


DH Tau B" 


JHK(1.10-2.46/mi) 


2600 


4.0 


2400 


3.0 2000 


4.0 


2350 


3.7 


>2400* 


4.0 


2 


J(1.10-1.34/im) 


2400 


4.0 


2200 


3.5 1600 


3.5 


2350 


4.0 


>2400* 


4.0 


2 


HK (1.46-2.46//m) 


2600 


3.5 


2500 


3.5 2700 


4.0 


2550 


3.8 


>2400* 


4.0 


2 


GQ Lup B" 


JHK(1.10-2.46jum) 


1900 


3.5 


1900 


3.5 1700 


3.5 


2050 


3.5 


1500 


4.0 


1 


J(1.10-1.34jum) 


2500 


4.0 


2400 


3.5 2600 


4.0 


2500 


4.0 


>2400* 


4.0 


2 


HK(1.46-2.46yum) 


1700 


4.0 


1700 


3.5 1700 


4.0 


1750 


4.0 


1500 


4.0 


1 


CT Cha B" 


JHK(1.10-2.46//m) 


1900 


4.0 


1600 


4.0 1700 


3.5 


1950 


4.0 


1500 


4.0 


1 


J(1.10-1.34yum) 


1900 


3.5 


1400 


4.0 1700 


3.0 


2000 


3.9 


1000 


4.0 


1 


HK (1.46-2.46yCim) 


2500 


3.5 


2300 


3.0 1800 


4.0 


2350 


3.1 


2400 


4.0 


2 


8 Myr Objects 


2MASS 1207 A 


JHK(1.10-2.46/mi) 


3200 


4.0 


2800 


4.5 >3000" 


4.5 


2900 


4.5 


>2400* 


4.5 


1 


J(1.10-1.34yum) 


2800 


4.5 


2700 


4.5 2900 


4.5 


2800 


4.5 


>2400* 


4.0 


2 


HK (1.46-2.46yum) 


3100 


4.5 


3000 


3.0 >3000 6 


4.5 


3050 


4.5 


>2400* 


4.5 


1 


TWA 5 B 


JHK(1.10-2.46/mi) 


2500 


3.5 


2300 


3.0 2200 


3.0 


2450 


3.3 


>2400* 


4.0 


2 


J(1.10-1.34^m) 


2700 


4.5 


2600 


4.0 2800 


4.5 


2700 


4.5 


>2400* 


4.0 


2 


HK (1.46-2.46ytim) 


2500 


3.5 


2000 


3.0 1700 


4.5 


2250 


3.1 


1600 


4.5 


1 


2MASS 1207 B 


JHK(1.10-2.46/mi) 


1600 


3.5 


1500 


3.5 1500 


4.0 


1650 


4.3 


1300 


4.0 


1 


J(1.10-1.34yum) 


1700 


4.5 


2000 


3.0 1700 


4.0 


1950 


4.0 


1500 


4.0 


1 


HK (1.46-2.46yum) 


1600 


4.5 


1500 


3.5 1500 


3.0 


1650 


4.5 


1300 


4.0 


1 


30-50 Myr Objects 


GSC 08047 B 


JHK(1.10-2.46/mi) 


2200 


3.5 
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Note: " These star formation region members have been dereddened as described in Sect. 4. 6 Note that this T cft falls at the edge of the grid used. 



Drift-PHOENIX, Gaia-Dusty, and Marley et al. The range in the 
parameters of T e g and log(g) are summarized in Table 5. The 
model fits were performed on each full J,H,K spectrum and the J 
and H+K subsets to investigate the variation of inferred physical 
properties on the wavelength coverage of the data. The results of 
the model fits for and log(g) are given in Table 6, divided by 
the age groups and ordered by average r eff from the J,H,K fit. 

Based on the J,H,K spectrum fits, the effective temperatures 
range from 1100 K to 3100 K and log(g) of 1.0 to 6.0, consid- 
ering the extremes of any given model fit. Some of the values 
for log(g) are outside the bounds expected from evolutionary 
models, and the fits also do not show a clear trend of increas- 
ing log(g) with increasing age. The lack of a clear correlation 
between log(g) and age may simply be a consequence of the 
coarseness of the grids and the dominant impact of temperature 
over surface gravity in shaping the overall spectral energy dis- 
tribution, as opposed to specific gravity-sensitive features. The 
range of masses included in the sample also complicates the 
comparison of log(g) and age. 



6.2.2. Restricted log(g) Model Grid 

Since the targets all have well-constrained young ages, the spec- 
tra were also compared with a restricted model grid limited 
to log(g) values appropriate for each age range, as detailed in 
Section 4 and based on the evolutionary models shown in Figure 
2 (Chabrier et al. 2000). Further support for using a restricted 
log(g) grid comes from measurements of the Na index defined 
in (Allers et al. 2007) and shown to be sensitive to surface grav- 
ity, but independent of spectral type. For the six targets for which 
the J-band data had sufficient sensitivity (2M1207A, TWA 5B, 
AB Pic B, GSC08047 B, and 2M0141), five had indices ranging 
from 0.93 to 1.08, within the range associated with giant stars 
and young brown dwarfs in star-forming regions (Allers et al. 
2007). Only 2M0141 has a higher value of 1.14, though this tar- 
get may be older than the rest of the sample. The inferred values 
of r ei f from the restricted grid model fits to the data are used 
to estimate the range of possible temperatures consistent with 
the observations. Uncertainties based on comparisons of the five 
different synthetic atmospheres should be more representative of 
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Fig. 6. The SINFONI spectra ( black lines) and best-fit model (colored lines) to the full spectrum is plotted for the star-forming region targets: DH 
Tau B (left), GQ Lup B (middle), and CT Cha B (right). The restricted log(g) model grids are used. 
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Fig. 7. The SINFONI spectra (black lines) and best-fit model (colored lines) to the full spectrum is plotted for the TW Hydra member targets: 
2M1207 A (left), TWA 5B (middle), and 2M1207 B (right). The restricted log(g) model grids are used. 



the true range of values than uncertainties estimated from a com- 
parison with a single model set. 

The SINFONI data and best-fit atmosphere model (restricted 
grid) from each of the five models is shown in Figures 6, 7, and 
8 for the three age groups. The region of the spectrum most diffi- 
cult to fit for the models is the peak and slope on the blue side of 
the H-band. The models often show a steeper slope and higher 



peak than the data. The peak of the H-band is also shaped by 
the dust cloud opacity, which is largely flat across the H-band. 
A high cloud opacity can flatten the spectrum over this wave- 
length range, and may explain the model discrepancy in the H- 
band shape. The portion of the data with the lowest signal-to- 
noise is the blue edge of the J-band, and improving this section 
of the data would enable more reliable fits, as some of the tar- 
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Fig. 8. The SINFONI spectra (black lines) and best-fit model (colored lines) to the full spectrum is plotted for the Tucana-Horologium and Field 
targets: GSC 08047 B (left), 2M0141 (middle), and AB Pic B (right). The restricted log(g) model grids are used. 
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Fig. 9. The best-fit effective temperature inferred from a single bandpass (left: J-band, middle: H-band, right: K-band) portion of the spectrum of 
each companion as a function of the effective temperature based on the full J, H, K spectrum. Results from the full grid are shown with open circles 
and from the restricted log(g) grid with crosses. All five models give systematically higher temperatures from fits using only the J-band spectrum. 



gets show a range of model best fits with significantly different 
slopes across the J-band. AB Pic B presents the most unusual 
shape, with a J-band notably lower flux than the H-band and no 
model is able to match the combination of flux levels and spec- 
tral shape as well as the other targets without adding other effects 
artificially. Unlike 2M1207 B, the inferred radius is not anoma- 
lously low. The inferred radii of all targets is discussed in section 
6.2.3. 

The results of the model fits to the full and partial spectra of 
each target are given in Table 7. Although some of the effective 
temperatures based on fits to subsections of the data are differ- 
ent for the full and restricted grids, the T e ff values for all targets 
using the full J,H,K spectrum differ by 100K or less with either 
the full grid or the restricted log(g) grid. The exception to this is 
2M1207B where the difference in effective temperature is 200K. 
The trend of higher inferred temperatures from J-band fits com- 



pared to H+K-band fits is seen in both the full grid fits and the 
restricted log(g) fits, as shown in Figure 9. The higher values of 
T e ff from the J-band fits compared to the H+K fits are most pro- 
nounced for the lower temperature objects below 2200K. There 
may be too much dust in the models in this temperature regime, 
resulting in underestimated J-band flux which a J-band fit com- 
pensates for with a bias towards higher T e ff values. 

Figure 10 shows the inferred Teff from each model as a func- 
tion of the average r e ff value. The calculation of the average 
excluded any values that corresponded to the limits of a model 
grid, and those limits are noted in Figure 10. The range in best- 
fit temperatures for a given object is +150-300K, considering 
all of the five different models. While no model always sug- 
gests the highest or lowest temperature, the DUSTY and Gaia- 
Dusty models give typically higher effective temperatures, and 
the Drift-PHOENIX and Marley et al. models are usually asso- 
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Fig. 10. The best-fit effective temperature of each object inferred from 
the restricted log(g) fits to each model (DUSTY - grey, BT-Settl - green, 
Drift-PHOENIX - red, Gaia-Dusty - purple, Marley et al. - blue) plotted 
against the average of all the plotted model fits, excluding values at the 
edge of each grid. Considering all models, each object has a range of 
±150-300 K, and the dispersion in temperature does not show a clear 
trend as a function of average effective temperature. 
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Fig. 11. The radius as a function of age for the targets (filled circles) and 
evolutionary tracks (Baraffe et al. 2003) for several substellar objects: 
5 M] up (pink line), 20 M Jup (black dashed line), and 30 M Jup (blue line). 



ciated with the cooler temperature fits for each target, which is 
consistent with the typically lower amount of dust compared to 
the DUSTY models. 
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6.2.3. Inferred Radii 

The fitting procedure to obtain the best match for effective tem- 
perature from each model grid includes adjusting the object ra- 
dius until the flux level corresponds to the synthetic spectrum. 
For the restricted log(g) fits reported in Table 7, each model 
r e ff and log(g) combination has an associated object radius. The 
range of inferred radii from the five separate models is given in 
Table 8; values based on best-fits at an extreme of a model grid 
are excluded. Although the radius is estimated from a restricted 
log(g) atmosphere model grid, the radius is not fixed by the value 
of log(g). 

As expected for young objects in the process of contract- 
ing, the radii are larger than the ~0.1 R Q predicted for older 
field planets and brown dwarfs (Baraffe et al. 2003). The only 
exception is 2M1207 B, and the apparent small radius - or un- 
derluminosity - of this object has been reassessed with updated 
models different from those considered in this study and found to 
be have an inferred radius consistent with evolutionary models 
(Barman et al. 2011). Due to the young ages of the targets, the 
radii may retain an imprint from the initial conditions, includ- 
ing the early accretion history of the objects (e.g., Baraffe et al. 
2002; Fortney et al. 2008; Baraffe & Chabrier 2010). Figure 11 
plots the inferred radii of the targets as a function of age, along 
with the evolutionary tracks for the size of 5-30 M Jup objects. 
The data exhibit a systematic decline with stellar age and are 
largely within the range predicted by the models (Baraffe et al. 
2003). 



6.2.4. Metallicity Grid 

For one of the model grids - the Drift-PHOENIX set (Witte et al. 
2009) - a range of non-solar metallicities are also included, and 
the targets were fit with the expanded library of synthetic spectra 
to investigate the impact of metallicity on the inferred physical 
properties of T e s and log(g). The results of the metallicity com- 
parison are given in Table 9, including fits to both the full spec- 
trum and the J- and H+K-bands. For approximately half the tar- 
gets, the best fit metallicity over the full spectrum was the solar 
value, supporting the use of solar metallicities for these young 
stars. For the remaining targets with best fit non-solar metallici- 
ties, the impact on the effective temperature and surface gravity 
estimates was minimal and resulted in a change of one step size 
in log(g) or T e g, except in one case. For DH Tau B, the difference 
in best fit effective temperatures is 700K. Overall, the ability to 
fit for metallicity does not typically change the inferred proper- 
ties by an amount larger than the uncertainty estimated from the 
range in values from the different models. 
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Table 9. Best Fitting Atmosphere Model Parameters, Metallicity 
Comparison 



Wavelength Range Drift Drift, [M/H] =0.0 

r eff (K) logg [M/H] T eff (K) logg 

DH Tau B" 

JHKQ.10-2.46/oti) 2700 33 ^06 2000 40~ 

J(1.10-1.34/oti) 1600 3.5 0.0 1600 3.5 

HKQ.46-2.46/ffli) 2700 3.0 -0.6 2700 4.0 

GQ Lup B° 

JHKQ.10-2.46/mi) 1700 33 O0 1700 33 

J(1.10-1.34//m) 2600 3.5 -0.6 2600 4.0 

HKQ.46-2.46/mi) 1600 4.0 -0.6 1700 4.0 
CT Cha B" 

JHKQ.10-2.46/ffli) 1600 330 ^06 1700 33 

J(1.10-1.34/oti) 1700 3.0 0.0 1700 3.0 

HKQ.46-2.46/mi) 1800 4.0 0.3 1800 4.0 
2MASS 1207 A 

JHKQ.10-2.46/ffli) 2900 43 ^06 3000 5 43T~ 

J(1.10-1.34/oti) 2900 4.5 -0.6 2900* 4.5 

HKQ.46-2.46/mi) 2900 4.5 -0.6 3000* 4.5 

TWA 5 B 

JHKQ.10-2.46/ffli) 2200 330 0.0 2200 330 

JQ.10-1.34/oti) 2700 4.0 -0.6 2800 4.5 

HKQ.46-2.46/mi) 1800 4.0 0.3 1700 4.5 

2MASS 1207 B 

JHKQ.10-2.46/ffli) 1700 40 03 1500 40 

JQ.10-1.34/oti) 1500 3.5 -0.6 1700 4.0 

HKQ.46-2.46/ffli) 1700 4.0 0.3 1500 3.0 

GSC 08047 B ~ 

JHKQ.10-2.46/ffli) 1900 4.5 0.0 1900 4.5 

JQ.10-1.34/oti) 1600 3.5 0.0 1600 3.5 

HKQ.46-2.46/mi) 1800 4.5 0.3 1700 4.0 
2MASS 0141 

JHKQ.10-2.46/ffli) 1700 40 O0 1700 40 

JQ.10-1.34/oti) 1600 3.5 0.0 1600 3.5 

HKQ.46-2.46/mi) 1800 5.0 0.3 1800 5.0 
AB Pic B 

JHKQ.10-2.46/ffli) 1600 43 O0 1600 43 

JQ.10-1.34/oti) 1800 4.5 0.3 1600 3.5 

HKQ.46-2.46/mi) 1800 5.0 0.0 1800 5.0 
Note: a These star formation region members have been dereddened as 

described in Sect. 4. * Note that this T e g falls at the extremes of the grid 

used. 



6.3. Comparison with Previous Observations 
6.3.1. Effective Temperatures 

The inferred effective temperatures for the majority of the sam- 
ple are consistent with previous estimates, however two of the 
targets - CT Cha B and AB Pic B - are notably different. In 
both cases, the range of best fit values for Tes is shifted to cooler 
temperatures, with the range of T e g 1500-1950 K for CT Cha 
B compared to 2350-2850 K from previous estimates (Schmidt 
et al. 2008) and 1400-1800 K for AB Pic B compared to 1700- 
2100 K from the literature (Bonnefoy et al. 2010). The implica- 
tions for the masses implied from evolutionary models (Chabrier 
et al. 2000) are shown in Figure 12 for CT Cha B and in Figure 
13 for AB Pic B. The lower effective temperatures are associated 
with lower masses, overlapping with the planetary mass regime. 

The origin of the discrepancy in inferred T e g for CT Cha B is 
clear - the different treatment of extinction. By treating Ay as a 
free parameter rather than using the value measured for the pri- 
mary, Schmidt et al. (2008) find a much higher extinction which 
results in a large change in slope across the J-K range. For AB 
Pic B, a previous model comparison (Bonnefoy et al. 2010) con- 
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Fig. 12. The mass as a function of r cff for substellar objects, according 
to evolutionary models of 1 Myr and 5 Myr objects (Chabrier et al. 
2000). The range of effective temperatures for CT Cha B from the best 
fit models is shown with solid lines and the range previously reported 
(Schmidt et al. 2008) is shown with dashed lines for comparison. The 
age of CT Cha B is expected to be between 1 Myr and 5 Myr. 



sidered each of the J, H and K bands separately, and so the dif- 
ferent results reported here arise from the systematically higher 
temperatures derived from fits restricted to the J-band, as shown 
in Figure 6. The T e ff values in Table 7 estimated from the J-band 
portion of the spectrum are comparable to the upper bound of 
the previously reported effective temperature. In both cases, ob- 
taining spectra over a larger wavelength range will enable more 
accurate estimates of T e ff and determine if CT Cha B and AB Pic 
B are analogues of 2M1207 B and represent wide orbit compan- 
ions with planetary masses, ideal cases for high signal-to-noise 
measurements of very low mass companions. 

6.3.2. Variability 

For two of the targets, DH Tau B and 2M0141, there are pre- 
viously reported J,H,K spectra taken at different epochs from 
the SINFONI data, and the combination of data sets are used to 
test for variability in the spectral shape, such as observed for the 
young low mass object TWA 30 in the TW Hydra Association 
(Looper et al. 2010). DH Tau B and 2M0141 are also young 
objects bracketing the age of TWA 30, making these targets in- 
teresting comparisons to TWA 30. The first epoch of DH Tau B 
spectra were obtained in 2003 (K-band) and 2004 (J+H-band) 
with CISCO at Subaru (Itoh et al. 2005), several years prior to 
the SINFONI observations in 2007 (Table 3). Over the 3-4 year 
time period between observations, the spectrum shape does not 
change. The first epoch of the 2M0141 observations were ob- 
tained in 2004 with Spex on the IRTF (Kirkpatrick et al. 2006), 
two years prior to the 2006 SINFONI observations (Table 3). 
Like DH Tau B, no variability in the spectral shape is seen. It 
is not possible to search for a grey brightening or dimming of 
the spectrum, such as that reported in Apai et al. (2012, submit- 
ted) because the flux calibration of the SINFONI data relies on 
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Fig. 13. The mass as a function of T eff for substellar objects, according 
to evolutionary models of 10 Myr and 50 Myr objects (Chabrier et al. 
2000). The range of effective temperatures for AB Pic B from the best 
fit models is shown with solid lines and the range previously reported 
(Bonnefoy et al. 2010) is shown with dashed lines for comparison. The 
age of AB Pic B is expected to be between 10 Myr and 50 Myr. 



the same photometry used to calibrate the first epoch spectra,. 
Additional coordinated photometric and spectroscopic measure- 
ments are required to determine the frequency, wavelength de- 
pendence and physical origin of variability in low mass object 
atmospheres. 



7. Summary and Conclusions 

With data from the SINFONI integral field unit, an empirical 
grid of J,H,K spectra for young ~l-50 Myr substellar objects has 
been constructed and compared with a suite of five theoretical at- 
mospheric model grids - DUSTY (Allard et al. 2001), BT-Settl 
(Allard et al. 2003, 2010), Drift-PHOENIX (Helling et al. 2008; 
Witte et al. 2009), Gaia-Dusty (Rice et al. 2010), and Marley 
et al. (Ackerman & Marley 2001; Marley et al. 2003; Saumon 
& Marley 2008; Freedman et al. 2008). The spectra should also 
serve as observational comparisons for updated models includ- 
ing more advanced treatment of physical processes. Since most 
of these objects are companions, there is additional information 
on the age and extinction of the substellar object to aid in the 
model fitting. All spectra show the distinct triangular shape over 
the H-band, a region of the spectrum which was found to be par- 
ticularly difficult for most models to match in detail. 

Based on the best-fit synthetic atmosphere spectrum for each 
target from each restricted log(g) model grid, the uncertainty in 
the temperature estimate for a young, substellar object is + 150- 
300 K from fits to the J,H,K spectrum. This is considerably 
larger than the uncertainty usually inferred from comparisons 
to a single model grid (~ 50-100 K, cf. Burrows et al. 2006). 
For all targets except 2M1207B, the inferred radius from the fit 
is greater than 0. 1 R , and generally consistent with the expec- 
tations of evolutionary models of young brown dwarfs (Baraffe 
etal. 2003). 



We note that fits to small wavelength ranges do not account 
for the broad spectral morphology and so may suggest effec- 
tive temperatures incompatible with observations of a broader 
range in wavelength. Fitting only the J-band portion results in 
systematically higher temperatures than fitting the H+K-band, 
likely due to systematic uncertainties in the model dust opac- 
ities. When the spectrum covers the full J,H,K range, the in- 
ferred temperatures are approximately the same whether or not 
the model grids include surface gravities limited to the range for 
young objects or expanded to include the higher values of older 
objects. While the temperatures were consistent using the full or 
restricted log(g) model grids, it is important to note that the spec- 
tral shape fit to the full model grid did not find the expected low 
surface gravities of these objects with known young ages. This 
apparent inconsistency is due to the fact that differences in spec- 
tral morphology caused by different effective temperatures dom- 
inate over spectral changes due to low surface gravity. Empirical 
comparison with observed spectra such as these may provide a 
more reliable method to determine the surface gravity; specific 
spectral lines and indices will be investigated in a follow-up pa- 
per (King et al. 2011). For the Drift-PHOENIX model grid, it 
was possible to include fits to non-solar metallicities, and vary- 
ing the metallicity did not impact the inferred temperatures by 
more than the uncertainty estimated from the range of model fits 
presented here. 

From comparison of the SINFONI data with previous spec- 
tra of two targets, there is no evidence of substantial changes in 
the spectral shape, however more observations and coordinated 
photometry and spectroscopy are required to rule out variabil- 
ity in these young substellar objects. The spectrum of AB Pic 
B shows an unusually low J-band spectrum relative to the H- 
band spectrum, analogous to 2M1207B, and the overall shape 
of the spectrum is more difficult to fit with the model grid than 
the remaining targets. For CT Cha B and AB Pic B, the new fits 
to the spectra suggest that the objects may have lower tempera- 
tures and, consequently, lower masses compared to previous es- 
timates (Schmidt et al. 2008; Bonnefoy et al. 2010). Additional 
wavelength coverage in the spectra will refine the temperature 
estimates to determine if these companions represent additional 
examples of planetary mass companions that can be investigated 
with high signal-to-noise spectroscopy, like 2M1207 B. 
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